function sim = set_up_mit_entry(eq_SS, param, glob, options, sim_in)

% entry shock
    sim = sim_in;
%     sim.D = sim.D/sim.D(end)*eq_SS.D;
% 
%     sim.A          = zeros(options.T, 1);
%  %   sim.A(1)       = -.01;
%  %   sim.A(2:100)   = sim.A(1)*(.5.^(2:100));
%     sim.A          = exp(sim.A);
%        
    rho = .685;
    
    sim.cE         = param.ce*ones(options.T,1);
    sim.cE(1:100)  = param.ce*(1+.001*(rho).^((1:100)-1));

    
%     sim.cF         = param.mu_f*ones(options.T,1);
% 
%     
       sim.C          = ones(options.T, 1);%.*((param.ce./sim.cE));
       sim.W          = sim.C;
%     
% %    sim.D          = eq_SS.D*ones(options.T,1);%.*((sim.cE./param.ce));
% 
%    % sim.D          = (sim.D-eq_SS.D)*.25+eq_SS.D;
%        
%     sim.mu         = nan(size(eq_SS.L,1),options.T);
%     sim.mu(:,1)    = eq_SS.L;
%     
%     sim.y          = nan(size(eq_SS.L,1),options.T);
%     sim.p          = nan(size(eq_SS.L,1),options.T);
%     sim.bp         = nan(size(eq_SS.L,1),options.T);
%     
%     sim.enter      = nan(size(glob.agrid,1), options.T);
%     sim.exit       = nan(size(eq_SS.L,1),options.T);
    
end